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Abstract. Modified Gravity (MOG) has been used successfully to explain the rotation curves of 
galaxies, the motion of galaxy clusters, the Bullet Cluster, and cosmological observations without 
the use of dark matter or Einstein's cosmological constant. We review the main theoretical ideas and 
applications of the theory to astrophysical and cosmological data. 
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1. INTRODUCTION 

The ingredients of the standard model of astrophysics and cosmology are: 

1. General Relativity, 

2. Large-scale homogeneity and isotropy, 

3. 5% ordinary matter (baryons and electrons), 

4. 25% dark matter, 

5. 70% dark energy, 

6. Uniform CMB radiation, T ~ 2.73 degrees, 

7. Scale-free adiabatic fluctuations AT /T ~ 10^^. 

Although the model fits available astrophysical and cosmological data [1], it opens up 
the mystery that about 95% of all matter and energy are invisible. The dark matter (DM) 
does not interact with ordinary baryonic matter and light. No current experiments have 
succeeded in detecting DM . The SNla supernovae data [2, 3] have created the need for 
the expansion of the universe to accelerate, promoting the mechanism of dark energy. 

In the event that DM is not detected, then to fit all available astrophysical and cosmo- 
logical data, we are required to modify Newtonian and Einstein gravity without assum- 
ing the undetected DM. A modified gravity (MOG) theory, also known as Scalar- Tensor- 
Vector Gravity or STVG [4, 5, 6], is based on an action that incorporates, in addition to 
the Einstein-Hilbert term and the matter action, a massive vector field, and three scalar 
fields corresponding to running values of the gravitational constant, the vector field cou- 
pling constant, and the vector field mass. 
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We begin in Section 2 by introducing the theory through the action principle, and 
establish key assumptions that allow us to analyze physically relevant scenarios. In Sec- 
tion 3, we derive the field equations using the variational principle. In Section 4 we solve 
the field equations in the static, spherically symmetric case. In Section 5, we postulate 
the action for a test particle, and obtain approximate solutions to the field equations for 
a spherically symmetric gravitational field. In Section 6 we demonstrate how the Fried- 
mann equations of cosmology can be obtained from the theory. In Section 7, we utilize 
the theory to obtain estimates for galaxy rotation curves, galaxy cluster dynamics and 
show how the solutions we obtained for the field equations remain valid from cosmo- 
logical to solar system scales. Lastly, we end in Section 8 with conclusions. 



2. MODIFIED GRAVITY THEORY 



The action of our theory is constructed as follows [4]. We start with the Einstein-Hilbert 
Lagrangian density that describes the geometry of spacetime: 



1 



{R + 2A),/^, 



(1) 



where G is the gravitational constant, g is the determinant of the metric tensor g^y (we 
are using the metric signature (+, — , — , — )), and A is the cosmological constant. We set 
the speed of light, c— I. The Ricci-tensor is defined as 



R 



(2) 



where F^y is the Christoffel-symbol, while R = g^^R^i^. 

We introduce a "fifth force" vector field via the Maxwell-Proca Lagrangian den- 
sity: 



-g, 



(3) 



where B^y = d^(j)v — dy^fi, jl is the mass of the vector field, co characterizes the strength 
of the coupling between the "fifth force" and matter, and is a self-interaction potential. 

Next, we promote the three constants of the theory, G, jU,(0, to scalar fields by 
introducing associated kinetic and potential terms in the Lagrangian density: 
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G 
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where denotes co variant differentiation with respect to the metric g^v, while Vg> Vju, 
and Vq) are the self-interaction potentials associated with the scalar fields. 
Our action integral takes the form 



(5) 
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where is the ordinary matter Lagrangian density, such that the energy-momentum 
tensor of matter takes the form: 



2 55m 



(6) 



where Sm = J d'^x. A "fifth force" matter current can be defined as: 



y 1 SSm 



r = 



(7) 



We assume that the variation of the matter action with respect to the scalar fields 
vanishes: 

0. (8) 



where Z = G,jU,(0. 



5X 



3. FIELD EQUATIONS 



The field equations of the theory can be obtained in the form of the first and second-order 
Euler-Lagrange equations [5]: 
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2^v J^M 



d(j)v 



J'' 



(9) 



VVVvM - - + ^Go^V?^^^^ + ^Vm(m) - = 0, (10) 



+v;(ft))=0, (11) 



G 



-y^(G) + G 



+ ^(i? + 2A) = 0, 



(12) 
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(13) 



4. STATIC, SPHERICALLY SYMMETRIC VACUUM SOLUTION 



In the static, spherically symmetric case with line element 



ds^ = Bdt^ -Adr^ - r^dO}, 



and with dQ? = dd^ + sin^ dd(j)'^, the field equations are written as 
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(18) 
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where 



// 2 , G' 1 ^ 2^x2 .2\ 1 fB' A'\ , 



+ 8^^'^'" " -^^«(«) = 0, (19) 
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inGT! = -A-V--N + ^^-^ + -, + - + -- 
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SnGTS = SnGT^ = -A-V - -N + 

^ ^ A 2A^r AA^B 

2AB?^4A^~2AB^AG^~i^AG^ AG^~2A^ 



IB" 


G" 
+ AG 


1 G' 

^ rAG 


2AB 


IB'G' 








-(ogI 


2ABG 


, AB 



B 



+ ^], (23) 



~AB~ 2AB^ ~ 2A^B ABr~ TJ^r^ Ap-~ ^ ^ 

^ = -4^ (^ + -^ (25) 



y = 2(oGV0((|)) + 8;r 



(26) 



G2 ' ;U2 

The prime denotes differentiation with respect to r, i.e., / = dy/dr. 

These equations can be substantially simplified in the matter vacuum case {Ty — 
0), with no cosmological constant (A = 0), setting the potentials to zero (V^ —Vg = 
V]u, V(o = 0) and also setting 0;- = 0. These choices leave us with six equations in the six 
unknown functions A, B, <j>t, G, [l, and ft): 

+ 2ft)G I^ + ^lJl- --— + ■ 



2ABG AGr \AB B J 2AB 4AB^ 

A'B' B' A' 1 1 



AA^B 2ABr 2A^r Ar^ 



+ ir^-^ = 0, (27) 
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G" 2G'2 B'G' A'G' B' A' 
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«G -7^ - +4;r + — 



A5 5 y \AG'^ Aju2 A 

B'G' 2G' B' 1 1 
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2A5G AGr ABr Ar^ 

,, 2 , u'2 G' 1 /5' A'\ , A(oG(b} , ^ 

// 2 , G' , 1/5' A'\ , G ^<^M^0f^ « 

G"+aG'4^+ir^-^^G'+ir^C-«'^^G 



r 2G 2V5 A ) 2VM^ 

--^AGi? = 0, (32) 

The values of A, 5, and B' are fixed by the requirement that at large distance from 
a source, we must be able to mimic the Schwarzschild solution (albeit with a modified 
gravitational constant), and that at spatial infinity, the metric must be asymptotically 
Minkowskian. The vector field must also vanish at infinity, which provides another 
boundary condition. Next, we assume that the values of G, fi, and (O are dependent on 
the source mass only, i.e., G' = fl' = (o' = 0. We seek the remaining initial conditions in 
the form of the fifth force charge Qs, and initial values of G = Go, /i = /io, and co = cOq. 
We note that the basic properties of the numerical solution and the solution's stability 
are not affected by the values chosen for these parameters. However, their values must 
be chosen such that they correctly reflect specific physical situations. To determine these 
values, we now turn to the case of the point test particle. 



5. TEST PARTICLE EQUATION OF MOTION 

We begin by defining a test particle via its Lagrangian: 

j^Tp = -m+ a(O^50juM^, (33) 

where m is the test particle mass, a is a factor representing the nonlinearity of the theory 
(to be determined later), 03 is present as it determines the interaction strength, ^5 is the 
test particle's fifth-force charge, and = dx^/ds is its four- velocity. 
We assume that the test particle charge is proportional to its mass: 

qs = Km, (34) 
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with K" constant and independent of m. This assumption implies that the fifth force charge 
^5 is not conserved, as mass is not conserved. This is the case in Maxwell-Proca theory, 
as ^^Jjx 7^ 0. We also have that the fifth force source charge Q=, <x M. 
From (33), the equation of motion is obtained 



m 



^^ + r^^M"M^^ = -a/cwmfi^vM"- (35) 



That m cancels out of this equation is nothing less than a manifestation of the equivalence 
principle. 

Our acceleration law can be written as [5]: 



r = 



GnM 



[l + a-a{l+iir)e-^"'], (36) 



where Gn is Newton's gravitational constant and a is given by 



a= J" (37) 



where £ is a constant of integration. 

The acceleration law (36) can also be recast in the commonly used Yukawa form: 



GyM 

—;t- 



l + ar(n-0e-'-/^ 



(38) 



with the Yukawa parameters Uy and X given by 



ay = iG~-G.)M 

{G^-GN)M + GN{y/M + E)^ 

X = l/M = ^. (41) 

Here, E and D are two universal constants of integration which can be determined from 
fits to galaxy rotation curve data. 

We can also express the acceleration law (36) as 

.= -^^, (42) 
where the effective gravitational constant Ggff is defined as 

Geff=Giv[l + a-a(l+jUr)e-^'']. (43) 
The metric parameter B{r) is given by 

= + (44) 

r 

The B{r) and A(r) solutions are shown in Figure 1. 
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FIGURE 1. Comparing MOG numerical solutions to the Reissner-Nordstrom solution, for a lO'' Mq 
source mass. The MOG metric parameters A (solid red line) and B (dashed brown line) are plotted along 
with the Reissner-Nordstrom values of A (dash-dot green line) and B (dotted blue line). Horizontal axis is 
in pc. We observe that the A metric parameter reaches at below the Schwarzschild radius of a 10^' Mq 
mass, which is ^ 0.01 pc. 



6. COSMOLOGY 



In the case of a homogeneous, isotropic cosmology, using the Friedmann-Lemaitre- 
Robertson-Walker (FLRW) line element. 



ds^ = dt^ -a^{t)[{l- k?)-^dr^ + ?dQL\ 
the field equations assume the following form: 

G 1 1 

CO + 3H(b - -G) - —G^^^l + — GV^ + = 0, 



(45) 

(46) 
(47) 



(0 



G ^ 3G / d , ^ 
+— A-— -+H^ 



0, 



(48) 
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(49) 
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where H = d/ais the Hubble expansion rate. 

It is possible to obtain an exact numerical solution to this set of equations using 
numerical methods [5, 6]. To carry out the solution, we assume a pressureless matter 
equation of state w = p/p = 0. Detailed fits to cosmological data, including the CMB 
angular power spectrum, the matter power spectrum and the SNla supernovae data have 
been obtained [6]. We find that the solutions can yield a "bouncing" cosmology. The 
bounce can be fine-tuned by choosing an appropriate value for Vg- This ensures that 
the universe reaches sufficient density in order to form a surface of last scattering. We 
emphasize that in our model only ordinary baryonic matter is present with a matter 
density of ~5% of the critical density. Nevertheless, the cosmology is flat, due in part 
to the increased value of the gravitational constant G, and in part to the presence of the 
non-zero energy density associated with Vq. 



7. FITTING GALAXY, CLUSTER DATA AND SOLAR SYSTEM 

DATA 

Unless one assumes that a massive dark matter halo is present, a typical spiral galaxy 
is dominated in mass by the central bulge. The motion of stars in the outer reaches of 
a galaxy can, therefore, be well approximated by the equations of motion in a static, 
spherically symmetric vacuum field. Indeed, our experience shows that the flat rotation 
curves of galaxies provide a sensitive test to determine the values of the constants D and 
E. In particular, it is easy to see that our results so far are compatible with the TuUy- 
Fisher law [10]. 

Kepler's laws of orbital motion yield a relationship between circular orbital velocity 
Vc at radius r from a mass M in the form 

!i = ^??. (52, 

TuUy and Fisher [10] have determined that for galaxies, assuming that the brightness of 
a galaxy and its mass are correlated, the flat part of the rotation curve obeys the empirical 
relationship: 

v^ocM, (53) 
where 3 < n < 4. In our case, we obtain 

V? °c Vm, (54) 
corresponding to n = 4 in the TuUy-Fisher relationship. 
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FIGURE 2. Photometric fits to galaxy rotation curves. There are 2 benchmark galaxies presented here. 
Each is a best fit via the single parameter (M/L) based on the photometric data of the gaseous (HI plus 
He) and luminous stellar disks. The radial coordinate (horizontal axis) is given in kpc and the rotational 
velocity (vertical axis) in km/s. The red points with error bars are the observations, the solid black line is 
the rotation curve determined from MOG, and the dash-dotted cyan line is the rotation curve determined 
from MOND [7]. The other curves are the Newtonian rotation curves of the various separate components: 
the long-dashed green line is the rotation curve of the gaseous disk (HI plus He) and the dotted magenta 
curve is that of the luminous stellar disk (from [8, 9]. 



Taking the next step, we select a small sample of galaxies and obtain an approximate 
fit to these galaxies yielding the values 

D ~ 625OM0/\pc"\ (55) 
E ~ 25000 Mg/l (56) 

The galaxy rotation curves we obtain for galaxies of varying mass are in good agree- 
ment with these values, treating D and E as universal constants without dark matter (Fig- 
ure 2). The galaxy rotation curves were obtained modeling the galaxies as point masses, 
benefiting from photometric data, as in the more extensive fit to galaxy rotation veloci- 
ties [8, 9, 11]. This exercise demonstrates that our established relationships between M, 
a, and jU not only satisfy the TuUy-Fisher relationship, but also offer good agreement 
with actual observations. N-body simulations of galaxy rotation curve dynamics using 
MOG have been performed [12]. 

In [9, 13], the spherically symmetric, static vacuum solution was used successfully 
to model galaxy clusters. We are able to produce a comparable result, while keeping 
the parameters D and E constant, by introducing an additional assumption: that the 
values of the MOG parameters Gcx, and ji at some distance r from the center of a 
spherically symmetric distribution of matter are determined not only by the amount of 
matter contained within radius r, but by the amount of matter within radius r*. Figure 3 
shows the case of r* = 3r. 

We have also succeeded in fitting the bullet cluster data [15], using MOG without dark 
matter [16, 17]. 
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FIGURE 3. A small sample of galaxy clusters studied in [9, 13]. Thin (black) solid line is the mass 
profile estimate from [14]. Thick (blue) solid line is the mass profile estimated using our STVG results. 
Dashed (green) line is the result published in [13], while the dotted (red) line is the Newtonian mass profile 
estimate. Radial distances are measured in kpc, masses in Mq. 



We have applied IVIOG to predict dispersion velocity curves for globular clusters, 
and found that the predictions follow those of Newtonian gravity [18]. By using Sloan 
Digital Sky Survey (SDSS) data, we have investigated how modified gravity theories 
including IVIOND and IVIOG affect satellite galaxies with the result that the data cannot 
currently differentiate significantly between modified gravity theories and dark matter 
models [19]. The IVIOG prediction for lensing caused by intermediate galaxies and 
clusters of galaxies has been investigated [20] . 

The theory must also be consistent with experiments performed within the solar sys- 
tem or in Earthbound laboratories. Several studies (see, e.g., [21]) have placed stringent 
limits on Yukawa-like modifications of gravity based on planetary observations, radar 
and laser ranging, and other gravity experiments. However, our prediction of the abso- 
lute value of the ay parameter is very small when Ay is small. The latter is estimated at 
Ay ~ 0. 16 pc (~ 5 X 10^5 m, or about 33,000 AU) for the Sun, and Ay ~ 2.8 x 10"^ pc 
(~ 8.7 X 10^^ m, or ~ 58 AU) for the Earth. The corresponding values of |ay| are 
|ay| ~ 3 X 10^^ and |ay| ~ 9 x 10^^^, respectively, clearly not in contradiction with 
even the most accurate experiments to date (Figure 4). 
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FIGURE 4. Predictions of the Yukawa-parameters from the MOG field equations are not in violation of 
solar system and laboratory constraints. Predicted values of A (horizontal axis, in m) vs. jay | are indicated 
by the solid red hne. Plot adapted from [21]. 



In the solar system the MOG field equations become essentially those of the Jordan- 
Brans-Dicke model [22, 23], for the influence of the vector field is reduced to very 
small values as shown in Figure 4. However, the standard JBD model coupling constant 
(DjBD has to be fine-tuned (Ojbd > 40, 000 to fit the Cassini spacecraft measurement of the 
Eddington-Robertson, parameterized post-Newtonian parameter 7 — 1 = (2.1 ±2.3) x 
10^; the other parameter /3 satisfies /3 = 1 in MOG. We have resolved this problem 
in MOG by coupling the scalar field G directly to matter by means of a scalar matter 
current: 

J = -\gT, (57) 

where T = T^^. This leads to obtaining an agreement with Earth based equivalence 
principle experiments and 7=1 [24] . 

We have plotted M vs. ro = jl^ in Figure 5. For the purposes of this plot, we 
used previously published results, while noting that our new calculations place dwarf 
galaxies, galaxies, and galaxy clusters by definition exactly on the line representing our 
prediction. This plot demonstrates the validity of MOG from the scales of star clusters 
to cosmological scales. 

We have investigated the possibility that MOG can explain in a fundamental way 
the origin of inertial mass. The static, spherically symmetric solution does not satisfy 
Birkhoff's theorem as in the case of the Schwarzschild solution in GR. This leads to a 
Mach-type influence of distant matter that can determine the inertial mass of a body. A 
possible spacecraft experiment has been proposed to test this prediction [25]. 

On the scale of Earth-based laboratory and solar system experiments with ever greater 
precision, MOG may eventually be verified or falsified. Beyond the solar system, as 
larger galaxy samples become available, the presence or absence of baryonic oscillations 
in the matter power spectrum may unambiguously decide in favor of modified gravity 
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FIGURE 5. The relationship jX^M = const, between mass M and the Yukawa-parameter ro — n ^ across 
many orders of magnitude remains valid. The solid red Une represents our theoretical prediction. 



theories or dark matter [5,6]. Confirmed detection of dark matter particles in deep space 
or in the laboratory would also be a strong indication against modified gravity. 



8. CONCLUSIONS 

In this paper, we have demonstrated how results of our Modified Gravity (MOG) theory 
can be derived directly from the action principle, without resorting to the use of fitted 
parameters. After we fix the values of some integration constants from observations, no 
free adjustable parameters remain, yet the theory remains consistent with observational 
data in the two cases that we examined: the vacuum solution of a static, spherically sym- 
metric gravitational field, and a cosmological solution. These solutions were explored 
using numerical methods, avoiding the necessity to drop terms or make other simplify- 
ing assumptions in order to obtain an analytic solution. Further, the constraints used to 
compute the solutions are consistent to the extent that they overlap with one another. 
The fact that at the level of the calculations presented here, our theory is not obviously 
falsified is an indication that we should pursue MOG further, for instance by obtaining 
interior solutions to the MOG field equations, and using these solutions to develop tools 
to perform A^-body simulations. 
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